Zero vorticity condition in calculation of ground state energy of Josephson junction 

lattices 

Y. Azizi 

Institute for Advanced Studies in Basic Sciences, Zanjan 45195-1 159, Iran 

We employ the charge neutrality condition in the form of zero vorticity for a Josephson junction 
lattice in calculating its ground state. We consider a Fibonacci ladder and the Penrose lattice to 
test our method. We also compare the results with those of the model that treats the plaquettes 
independently. We show that the zero vorticity condition improves on the results of the independent 
plaquette model. 
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I. INTRODUCTION 

A Josephson junction is characterized by its current- 
phase relation I((j>), where <f> is the phase difference across 
the junction. I(4>) must have two general properties: it 
must be 2ir- periodic and it must be an odd function. 
Therefore I(cf>) can be any periodic, odd function^— . 
In here we focus on sinusoidal Josephson junctions. In 
this case, the Hamiltonian of system is a cosine func- 
tion of the gauge invariant phase difference^. A Joseph- 
son junction lattice is a lattice with nodes that are su- 
perconducting islands and bonds that represent Joseph- 
son junctions. This lattice of Josephson junctions is a 
model that can be used for the description of other im- 
portant systems like stacked Josephson junctions^, high- 
temperature superconductors^—, and a rotating BEC 9 . 

If we ignore screening effects, we can write the Hamil- 
tonian of the system as a sum of Hamiltonians of Joseph- 
son junctions^. In here we focus on minimization of this 
Hamiltonian as a function of superconducting phase in 
presence of external magnetic fielcU^. By this minimiza- 
tion we obtain the ground state energy as a function of 
the external magnetic field, or frustration factor, which is 
defined as the ratio of flux through a plaquette with unit 
area to the quantum of flux. Behavior of ground state 
energy E as a function of frustration factor / is interest- 
ing. When / is zero all junctions can attain their abso- 
lute minimum, but when / becomes nonzero, the bonds 
of the lattice cannot generally obtain this absolute mini- 
mum simultaneously and produce some superconducting 
current-this is the origin of the term frustration. 

We approximate the behavior of the lattice in ground 
state by the Independent Plaquette Model (IPM)^ 
which considers the plaquettes independently. Also we 
suppose that total vorticity of lattice is zero. This con- 
dition can be understood once we look at the relation 
of this model with Charged Coulomb Gas (CCG) model. 
The Hamiltonian of the Josephson junction lattice can be 
mapped to CCG on dual lattice, where the charge of each 
point in dual lattice is proportional to vorticity of cor- 
responding plaquette in original lattice^— . This shows 
that the zero vorticity condition in CCG model means 
that the system is neutral. 

This condition was applied before for the description 



of numerical results on the square lattic e) 11 ' 12 . In here 
we apply this condition for lattices with two types of 
plaquettes, corresponding to CCG models with four types 
of charge. Our result is in agreement with those without 
this explicit condition and is an improvement on thenJ^. 

Behavior of E(f) is in close relation with behavior of 
critical temperature Tc(f) of lattice. In fact, when the 
ground state energy is increased, we expect that the crit- 
ical temperature is lowered and vice versa^i . Our result 
is in agreement with the mean field approach for Tc(f)^. 

This model can be mapped onto the frustrated XY 
mode l 10 ' 15 ' 17 , which is the XY model with unitary cou- 
pling which is dependent on a parameter equivalent to 
the frustration factor. This model is in close relation to 
the Frenkel-Kontorova mode l 18 ' 19 which describes the be- 
havior of the ground state for a system of particles under 
the effect of the lattice potential. 

The structure of the paper is as follows: in section II, 
we introduce our model with the zero vorticity condi- 
tion; in section III we compare the results of the model 
with the numerical study, comparing energies and corre- 
sponding Fourier spectra; the final section is devoted to 
Conclusions. 



II. THE MODEL 

For a Josephson junction lattice, the Hamiltonian is, 
H = -J2 cos (7«) (2-1) 

<i,j> 

where < i,j > means that ith node and jth node are 
connected through a bond of lattice. %j is the gauge 
invariant phase difference. It is equal to 

= Oi - Oj - A *,j ( 2 -2) 

where 9 is the superconducting phase and A^j is the in- 
tegral of the vector potential from ith node to jth node. 
We define the frustration factor as the ratio of magnetic 
flux through a plaquette with unit area to the flux quan- 
tum. 
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The gauge invariant phase difference satisfies the flux- 
oid quantization^, 

V k = Y J li3=M^k-a k f) (2.3) 

where n k is an integer and a k is the area of fcth plaque- 
tte. Summation is over the edges of fcth plaquette. T4 
is vorticity of plaquette, and in transformation to CCG, 
is proportional to charges on fcth site of the dual lattice 
pointii - — . 7ifc can have two values: |_ a fc/J or L a fc/J + 1 
where [^J means greatest integer less than x. It means 
that if we have I types of plaquettes, then we have 21 
different values for V% and therefore 21 types of charges 
in the equivalent CCG model. We want to minimize this 
Hamiltonian with respect to the superconducting phases. 
We use the Monte-Carlo (MC) method22i2I for the nu- 
merical minimization. We denote this minimum by E 
which is normalized by the number of junctions. 

From symmetry of a single plaquette, wc can say that 
all four gauge invariant phase differences are equal, and 
we denote it by 7. Then from Eq. 12. 31 we can find 
7 = 7r/2(n — af), for two choices of n, and giving the 
ground state energy of the plaquette as^S, 

E, = -cos(w/2([af\ + i - af)) (2.4) 

for i = 0,1. Now, we approximate the Hamiltonian 
as follows: supposing that fcth plaquette has vorticity 
[o-kf\ + ik, then it has energy Ef in IPM approxima- 
tion, and we can write Hamiltonian 12.11 as sum of these 
energies in IPM approximation: 

H IPM = Y,El (2.5) 

k 

where ik determines vorticity of fcth plaquette. Now we 
minimize above Hamilitonian with respect to i k - In orig- 
inal IPM, we choose ik, such that |_ a fc/J + becomes 
nearest integer to a,kf—- Because of interaction between 
plaquettes this is not exact, and some of the plaquettes 
give the furthest integer as their vorticity. We suppose 
that this deviation from IPM is such that total vorticity 
of lattice becomes as small as possible. This is the zero 
vorticity condition, with the total vorticity given as, 

Vtotai = ^ Vfc = 2n(n tota i - ai attice j) (2.6) 

k 

where n tot ai =J2k n k an d aiattice is the area of the whole 
lattice. The condition for zero total vorticity reads, 

Vtotal/N = J2 V */ N ~ ( 2J ) 

k 

where N is total number of plaquettes. 

This condition can be understood if we look at the rela- 
tion between Josephson junction lattice and CCG. It can 
be shown that Josephson junction hamiltonian can be 
mapped to CCG on the dual lattice, and in this trans- 
formation Vfc is proportional to charge on the equiva- 
lent point in the dual lattice. Hence, the zero vorticity 



condition is equivalent to the charge neutrality of CCG 
system»iir— 

Suppose that we have k types of plaquettes in a lattice, 
then each type can take on two energies, either Eq or 
Ei, bringing in 2k variables. But the total number of 
each type of plaquette is determined and therefore we 
have k constraints on these variables and there are k 
independent variables. Using the above equation lets us 
remove one of these variables and finally we have k — 
1 variables. If we apply the IPM approximation, the 
Hamiltonian becomes linear in terms of these variables 
and easily can be minimized. We denote this model as 
IPMO. 

We use the power spectrum of E(f) to characterize it, 

S(w) = EE* (2.8) 

where E* is the complex conjugate of E, the Fourier 
transform of E(f), 

E(u) = ±- [ E(f)exp(2iriLuf)df (2.9) 

The power spectrum was calculated in MATLAB by the 
FFT method^. 

Now we focus on lattices with two types of plaque- 
ttes. We choose two lattices: the Fibonacci ladder and 
the Penrose lattice. These two cases are in close relation 
with each other. Ratio of areas of plaquettes in both 
are the same and both make use of the Fibonacci order. 
Therefore we can compare their energies to give some in- 
sight about the difference between them. From the above 
discussion, we can find that these Josephson junction lat- 
tices are equivalent to the CCG model with four types 
of charge. 

III. LATTICE WITH TWO TYPES OF 
PLAQUETTES 

In here we consider the above formulation for lattices 
with two types of plaquettes which we denote them by 
L and S. S plaquette has area 1 and L plaquette has 
area a. Density of S plaquettes is N$ and that of the L 
plaquettes is Nl, therefore Nl + N$ = 1. Wc suppose 
that the density of S (L) plaquettes with energy Eq is 
ns(riL), therefore from Eq. 12.71 we can find that 

n s + n L = l + N L ([af\-af)+N s ([f\- f) (3.1) 

This condition for a lattice with one type of plaquette 
(aL = as = 1, ns = til = n) is reduced to 

»=1 + (L/J"/) (3-2) 

which is in agreement with the previous resul t 11 ' 12 . The 
Hamiltonian for a lattice with two types of plaquettes 
within the IPM is 

Hipm = {n S El + (N s - n s )E? + n L E^ + {N L - n L )E^) 

(3.3) 
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Now we minimize Hamiltonian with respect to these twc 
densities with condition Eq. 13.11 If we increase ul 
by one (therefore decrease ul by one), then change in 
Hamiltonian is AE = Efr - E{ - Eg + Ef , therefore ii 
AE > then n s = min(N s ,q(f)) and if AE < then 
iil = min(NL,q(f)) where min(x,y) means the mini- 
mum between x and y . We can write these results as 




n s = N s d(AE) + (q(f) - N L )(1 - 9(AE)) (3.4) 
n L = N L (1- 6{AE)) + (q(f) - N s )e(AE) (3.5) 



FIG. 2: Numerical and IPMO Fourier spectrum for the Fi- 
bonacci ladder 



where 9{x) is the step function which is zero for negative 
x and one for positive x, and q{f) represents the right 
hand side of Eq. 13.11 

Now we apply these results to a ladder with two types 
of plaquettes and the 2D Penrose lattice. 



A. Fibonacci Ladder 

For this lattice we have, a = r = (1 + \/(5))/2 and 
the plaquettes appear in the Fibonacci order, defined as 
follows. Denoting the structure in its mth step of con- 
struction by U m , the structure is constructed recursively 
by the rule U m+ 2 — U m+ i + U m , where summation means 
adjacent placement: to get the structure in step (to + 2). 
put the structure in step (to + 1), first and that of step 
to, to its right. This generates the sequence 

S, L, LS, LSL, LSLLS, LSLLSLSL, ... (3.6) 

In the limit of an infinite structure this defines a qua- 
sicrystalline lattice such that ratio of Nl to Ns is r. The 
energy and Fourier spectra for MC and IPMO are given 
in Figs. [I] and H 




FIG. 1: Numerical and IPMO energy for the Fibonacci laddei 

As we see, the agreement between two results is very 
good. From Fourier spectrum comparison, we see that 
there are peaks on u> = 1, r, 1 + r, r — 1, for both models. 
This agrees with the results of IPM without the zero 
vorticity condition^. 

We also compare results of these two methods for other 
types of lattices with different ratio of areas and again 
find good agreement. 



B. Two dimensional example: The Penrose lattice 

Penrose lattice is an example of 2D quasicrystals^ 
which has two types of rhombi with areas sin(2n/5) and 
sin(7r/5). It is known that minimization of Josephson 
Junction Hamiltonian on this lattice, as well as on the 
Fibonacci ladder, does not result in a periodic E(f); 
similarly, Tc(f) is aperiodio^i— . Our numerical results 
agree with this observation. 

We compare the ground state energy for this lattice 
with its one-dimensional counterpart, the Fibonacci lad- 
der. This shown in Fig. [31 We see that the ground state 
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FIG. 3: Comparison between energy for Penrose lattice Ep 
and Fibonacci ladder Ef- 

energy for the Penrose lattice is greater than that of the 
Fibonacci ladder. This can be related to the constraints 
involved: in the two dimensional case there are more 
neighboring plaquettes than in the ladder geometry. 

The ground state energy using MC and IPMO and their 
corresponding Fourier spectra for the Penrose lattice are 
shown in Figs. SjandO 




FIG. 4: Numerical and IPMO energy for Penrose lattice 
We can see that the results of IPMO are in agreement 

generally with Monte-Carlo results. Also IPMO predicts 

the main frequencies of E(f) correctly. 
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FIG. 5: Numerical and IPMO Fourier spectrum for Penrose 
lattice 

IV. CONCLUSION 



sum and difference. In 2D the accuracy of the model is 
not as good as in the case of the ladder, because of greater 
effects of neighboring plaquettes: In 2D a plaquette has 
more neighbors than in the ladder, and therefore IPM is 
worse in 2D than in the effectively ID ladder. 

Zero vorticity condition can be understood when we 
look at the connection of Josephson junction Hamiltonian 
with CCG. In CCG model vorticity of plaquette plays 
the role of charge on its corresponding dual lattice site. 
Therefore zero vorticity condition in CCG means that 
system of charge is neutral, which naturally leads to a 
lower energy for the ground state. 



We enter a condition in the calculation of the ground 
state energy of Josephson junction lattices, namely that 
the total vorticity of lattice be zero. This condition 
along with the Independent Plaquette Model (IPM) gives 
us the properties of the ground state energy with good 
approximation. It can describe the main peaks of the 
Fourier spectrum of E{f) on areas of plaquettes and their 
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